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Tb ■ Abstract 

^ , We study the ground state of a spin-half Heisenberg antiferromagnet on the stacked kagome 

lattice by using a spin-rotation-invariant Green's-function method. Since the pure two-dimensional 
(N' 

K^ ■ kagome antiferromagnet is most likely a magnetically disordered quantum spin liquid, we investigate 

CN , 
^^ . the question whether the coupling of kagome layers in a stacked three-dimensional system may lead 

js^ I to a magnetically ordered ground state. We present spin-spin correlation functions and correlation 

^— ^ ' lengths. For comparison we apply also linear spin wave theory. Our results provide strong evidence 

j^ ' that the system remains short-range ordered independent of the sign and the strength of the 

interlayer coupling. 

a 

O . PACS numbers: 75.10.Jni; 75.45. +j; 75.50.Ee 

> 

X 



I. INTRODUCTION 

In frustrated quantum spin lattices the interplay of quantum and frustration effects 
causes interesting physics, see e.g. the recent reviews |l||2|3|y|. Special attention was fo- 
cused on the problem of the ground state (GS) nature of the Heisenberg antiferromagnet 
(HAFM) on the two-dimensional (2D) kagome lattice. In the classical limit the GS of 
the HAFM on the kagome lattice exhibits a huge non-trivial degeneracy (see, e.g. Refs. 
|2|5||6|l7|) . Intensive work over the last decade on the spin half quantum version of the 
J^odel5,8,9,io,ii,i2,i3,i4,i5,i6,i7,i8,i9,2o,2i,22,23,24,25,26 jg^ ^^ ^j^g couclusiou that for the quautum 

HAFM the GS is most likely a short-range ordered spin liquid without any kind of long- 
range order (LRO). However, we mention that the GS properties in the quantum case are far 
from being fully understood, yet. Most of the conclusions are drawn from finite-lattice results 
of up to A^ = 36 sites. Very recently, the rotation-invariant Green's-function method intro- 
duced by Kondo and Yamaji^^ has been successfully applied on the spin-half HAFM on the 
2D kagome lattic o^^i^^ . In the framework of this approach one obtains also a spin-liquid GS 
with an energy per spin Eq/JN = —0.4296 which agrees well with the best available results 
obtained by other methods like the exact diagonalizatior4»2»i2*i^ {Eq/JN = —0.4344) and 
the coupled-cluster method^'^ {Eq/JN = —0.4252). Note that additional second-neighbor 
couplings may lead to a magnetically ordered GS phas o^^i^^i^^ . 

It is well-known that the dimensionality is crucial for the existence of order. Roughly 
spoken one can say, the higher the dimensionality of the spin system the more the influence 
of quantum or thermal fluctuations is suppressed. For the GS of quantum spin systems one 
finds several examples, where the physical properties change basically going from one to two 
dimensions or from two to three dimensions. For instance, the transition from the linear- 



chain to the square-lattice HAFM was studied in Refs. l30lJ3llJ32ll33l Another example is the 
frustrated J1-J2 HAFM, which exhibits a magnetically disordered quantum para magne tic 
GS phase around J2 ~ 0.5Ji for the 2D square lattice (see, e.g. Refs. yiaJadkUykd 
and references therein) but does not show this kind of quantum phase for the 3D body- 
centered cubic lattice^i^. Moreover, we know from the Mermin- Wagner theorem^ that the 
interlayer coupling in quasi-2D Heisenberg magnets is crucial for magnetic ordering at finite 
temperatures. The role of the interlayer coupling in quasi-2D systems becomes particularly 
interesting if the decoupled layers themselves are magnetically disordered, e.g. due to strong 





FIG. 1: Sketch of one kagome layer in the congruently stacked lattice with its in-plane geometrical 
unit vectors ai = (0, 2, 0) , a2 = (\/3, 1, O) . The out-of plane unit vector a.^ = (0, 0, 1) is not shown. 
Within a geometrical unit cell the spins are distinguished by a running index a = 1,2, 3. 

frustration. Though this question is relevant for real 3D solids with a layered arrangement 
of magnetic atoms it has been less studied yet. 

In this paper we address this question for the HAFM on a congruently stacked kagome 
lattice. For such a system it is not known, whether an appropriate interlayer coupling 
may lead to a semi-classically ordered magnetic GS. However, one can expect that in the 
special limit of infinite ferromagnetic interlayer coupling the physics of the classical 0(3) 
kagome antiferromagnet is obtained (see Sec. IV Bj) . We emphasize, that the Lanczos exact 
diagonalization of finite lattices being the most powerful method to study the HAFM in the 
pure 2D kagome case is inappropriate in 3D. Therefore we will discuss this question in our 
paper by the above mentioned spin-rotation-invariant Green's-function method, having in 
mind that this method was successfully applied to the pure 2D case. In particular, we will 
calculate the correlation functions and the correlation lengths. In addition, we will discuss 
the linear spin-wave theory (LSWT) for comparison. 

II. THE MODEL 



We consider congruently stacked kagome layers shown in FigC] The vertices of this 
stacked lattice are occupied by A^ spins one-half interacting via the Heisenberg exchange 



coupling 

-n ^ — y ^ Jma,nf3^ma^nf3j {^) 

ma,nf3 

where the sum runs over all unit cells (labeled by m and n) and all spins within a unit 
cell (labeled by running indices a and /3, see Fig. (T]). One unit cell contains three spins, 
therefore the number of cells is A^/3. The exchange coupling Jma,ni3 is non-zero for nearest 
neighbors (NN), only. We introduce two different exchange parameters Jy and J_l according 
to Jma,ni3 —* J\\ (1 — ^a,i3) + J±Sa,f3- Siucc wc are interested in the stacked kagome HAFM, 
we consider J\\ > but allow both signs for J±. 

III. CLASSICAL GROUND STATE AND LINEAR SPIN WAVE THEORY 

We start with a brief discussion of the classical ground state (GS) of ((H). Because there 
is no additional frustration caused by J^, the classical GS within a certain layer is not 
modified, and the nontrivial huge degeneracy due to corner-sharing triangles^j^*^ is not 
lifted. The only effect of a ferromagnetic (antiferromagnetic) interlayer coupling J± is the 
parallel (antiparallel) orientation of interacting spins of adjacent layers. Thus for a certain 
classical GS the spin configurations in all layers are identical. Two particular variants for 
the classical GS, namely the so-called -\/3x^/3 state and the so-called g = state (see e.g. 
Fig. 17 in Ref. |j|), are often used for discussing possible magnetic LRO in the kagome lattice. 
These two particular planar states can also be considered as variants of possible GS ordering 
for the stacked kagome lattice. 

Let us remind the reader of the fact, that starting from the v^x-\/3 state as well as from 
the g = state the linear spin-wave theory (LSWT) for the pure 2D kagome HAFM^*i^ 
leads to one fiat zero mode and two degenerate modes producing divergent integrals in the 
sublattice magnetization, which might be taken as a hint on the absence of semi-classical 
magnetic order. 

For the stacked (3D) system we start also from both the -\/3xa/3 state and the g = 
state and perform the LSWT as usual. Taking into account that we have three spins per 
unit cell we have to introduce three different kinds of magnons. As in the 2D case^^ the 
spin wave spectrum is equivalent for all coplanar classical GS configurations because the 
directional cosine of in-plane neighbors is always —1/2 while the directional cosine of out- 
of-plane neighbors is either 1 or —1 depending on the sign of J±. For J^ > we obtain for 



the spin-wave dispersions 



with 



^iq = s ^J AJl {1 - cos^q^) + 6 J\\J±{1 + COS Qz), (2) 

^2q = S^/a^T^, UJ3q, = S^/Oq - 6q (3) 

4 Jf (2 cos^ I + ( cos^ y^ + cos^ I ) (1 - 2 cos2 |^ 

+ 4Ji (1 - cos^ q,) + 3J\\Ji_ (3 - cosg^) , (4) 

K = J||J±(1 -3cosg^)Dq, (5) 

Dq = Jg + 16 cos4 ^ + 16 cos2 ^ cos2 ^ - 8 cos2 ^ - 24 cos^ ^. (6) 

^ V 2 2 2 2 2 ^ ^ 



Oq 



For J I < we find 



^iq = si74J|(l -cosg^) -6J||J_l(1 -cosg^), (7) 



^2q = SA/flq + feq, C^Sq = S A/Clq - &q (8) 



with 



aq = 4 Jf (2 cos2 I + ( cos2 ^ + cos2 | J (l - 2 cos^ |^ 

+ AJl (1 - cosg^)^ - 9J|| Jx (1 - cosg^) , (9) 

&q = J||J±(l-cosg,)Dq. (10) 

The quantity s in Eqs. Q , © > © > (0) is the spin quantum number which can be considered 
as parameter of the model © within the LSWT. If J± goes to zero, the known fiat zero 
mode as well as the both degenerate modes^*i^ are recovered. But even for non-zero J± we 
have a fiat zero mode c^iq and degenerate u;2q and c^sq in the q^ — qy plane for certain q^. 
For J_L > 0, ciJiq is a fiat zero mode at g^ = vr, and u;2q and cusq become degenerate at any 
qx, qy for g^ = arccos 1/3. For Jj_ < 0, c^iq becomes a fiat zero mode in the qx — qy plane for 
g^ = 0, and to'2q and tusq are degenerate at any g^, qy also for g^ = 0. Thus, in both cases 
there is no finite sublattice magnetization due to the resulting logarithmic divergencies in 
the involved integrals. We conclude, that the LSWT yields very similar behavior for the 
pure 2D and the 3D stacked case. 



For completeness we give the expression for the GS energy per spin Eq/N: 

3 

§ = - (j|| + 1 jj) 3 (3 + 1) + ^ E E ^-^' (11) 

q m=l 

where the Wmq are the respective spin wave dispersions from Eqs. ((21), © or ((Tj), (jH)). 



IV. SPIN-ROTATION-INVARIANT GREEN'S-FUNCTION THEORY 

The Green's function method is one of the most powerful techniques for the investigation 
of quantum many-body systems and was successfully applied to spin systems over many 
decades. The rotation-invariant decoupling scheme was introduced by Kondo and Yamaji^i 
to study the one- dimensional Heisenberg antiferromagnet and ferromagnet at finite tem- 
peratures. This decoupling scheme was developed to improve the description of magnetic 
short-range order and allows to describe magnetic order-disorder transitions driven by quan- 
tum fluctuations as well as by thermal fluctuations. Later on the spin-rotation-invariant 
Green's-function theory was used to discuss several two-dimensional models, for example 
the pure HAFM^^»^, the frustrated J1-J2 HAFM224& and the spatially anisotropic HAFM 
on the square lattice^ as well as the doped square-lattice HAFM (t-J-model)^^^. The 
quasi-two-dimensional and the three-dimensional HAFM have been investigated, to o^^i^° . 

Let us outline the main ideas of this method based on the equation-of-motion tech- 
nique. To evaluate the relevant correlation functions we have to calculate a set of Fourier- 
transformed Green's functions ((S'^Jq,; 5*1^^)) which are related to the dynamic spin sus- 
ceptibilities by Xqa/3 (^) = ~ (('S'qa; 5*1^)) . Their equation of motion reads 

^ ((^q".; SZ^))^ = ( K„ S-_^,]_) + ((z^i; SZ^,))^ . (12) 

Supposing rotational symmetry we have (S^^) = for any spin and, as a result, 
\ ['S'+„''S'Iq^] ) = 0. Furthermore, we have {S^^S^^) = \{S^^S~p), i.e. it is sufficient 
to calculate {S^^S~p) to know all components of the correlation function (SmaS^/j). Go- 
ing beyond the RPA decouplings^ we consider in a second step the equation of motion for 



uj{{iS^^;S_^p 



q+ C- 

^qa! "-"-q/g 



+ {{-SL:SZ,0)) . (13) 



The combination of Eqs. (fT^ and (fT^ yields 

{{sU szj)^ = ( [s^, 5:^J _) + ((-5i; 5:^,))^ . (14) 
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The operator —S^^ = [[S^^, H], H] contains products of three spin operators along NN 
sequences. Those operator products were treated in the spirit of the decouphng scheme 
by Shimahara and Takadai^. This decoupling is performed in the site representation 
of the Green's functions. For example, the operator product S^S'^S^ is replaced by 
Va,b (SaSb) Sc + Va,c (SaSc) S^, where A, B, C represent spin sites. The quantities 77^,^ 
are vertex parameters introduced to improve the approximation scheme. In the minimal 
version of the theory we introduce just as many vertex parameters as independent condi- 
tions for them can be formulated. Because there are just two such conditions (see below) we 
consider two different parameters r/y and rij_ attached to intralayer and interlayer correlators, 
respectively. After the decoupling we write Eq. (|T^ in a compact matrix form omitting the 
running indices a and /?, 

{u' - Fq) xt (^) = -^q' (15) 

where Fq and Mq are the frequency and momentum matrices, respectively. Since the unit 
cell contains three spins, F^, Mq and xt~ ^^^ 3 x 3-matrices. For the sake of brevity we 
do not give the lengthy expressions for the matrix elements of Fq and Mq but give their 
eigenvalues, only. Both matrices are hermitean and commute with each other. Hence, the 
solution of p5j) in terms of the common set of normalized eigenvectors |jq) of xt~ reads 

xt (^) = - E ^P^ i^'q) O'^i • (16) 



For the eigenvalues m^q of Mq and tu^q of Fq we find 

miq = -12J||Ci,o,o - 4J±Co,o,i (1 - cosg^) , (17) 

^iq = 3J|| (1 + 2?7|| (2ci,o,o + Ci,i,o + C2,o,o)) 

+Jl (1 - cosg^) (1 + 2?7_L (co,o,i + £0,0,2) - 4?7_lCo,o,i (1 + cosg^)) 

+4J|| J_L (r7±co,o,i (1 - cosg^) - 2?7_lCi,o,i (1 + cosg^) 

+ 7r7_LCi,o,i - 3r7||Ci,o,o cosg^) , (18) 



"^2q = -6J||Ci,o,o -4J±co,o,i (1 -cosg^) -2J||Ci,o,o-Dq, (19) 

'3 



^2q = 4 9 + 2^11 (2Cl,0,0 + Cl,l,0 + C2,0,o) 



I o ( o 4 ?y I o 2 V jo's; 2 Q'y 2 V -jo's „ 2 Q's/ 

+8r7||Ci,o,o 2 COS -^ + 2 cos — — cos -^ - cos — 3 cos -^ 

+Jl (1 - cosg^) (1 + 2r]± (co,o,i + 00,0,2) - 4r7_LCo,o,i (1 + cosg^)) 
+J\\J± (-877_lCi,o,i (1 + cosg^) - 2r]±cofl,i (1 - cosg^) 



6r7||Ci,o,o cos qz + 22r7_LCi,o,i) 
^ ( J|l (1 + 2r/|| (2ci,o,o + ci,i 
+ 4Jj_ (r7±co,o,i(l-cosg2) + r7_LCi,o,i - r7||Ci,o,o cosg^)) L'q, (20) 



+y (^11 (1 + 2^11 (2ci,o,o + Cl,l,o + C2,0,o)) 



"^3q = -6J||Ci,o,0 -4Jj_Co,0,l (1 -COSg^) +2J||Ci,o,0-Dq, (21) 

^3q = 4(2"^ ^^11 (2ci,o,o + Cl,l,o + C2,o,o) 

+8r;i|Ci,o,o 2 cos -^ + 2 cos cos -^ — cos 3 cos -^ 

+J]_ (1 - cosg^) (1 + 2r7_L (co,o,i + £0,0,2) - 4?7_lCo,o,i (1 + cosg^)) 

+J\\Jl (-877_lCi,o,i (1 + cosg^) - 2r7_LCo,o,i (1 - cosg^) 

- 6r7||Ci,o,ocosg2 + 22r7_LCi,o,i) 

~y (^11 (1 + 2^711 (2ci,o,o + Ci,i,o + C2,o,o)) 

+ 4Jj_ (r7±co,o,i(l-cosg2) + r7_LCi,o,i - ^||Ci,o,o cosg^)) ^q, (22) 

where Dq is defined in Eq. ^. The correlators Q,fc,m have to be determined self-consistently; 
their indices correspond to a vector R = /ai/2 + ka2/2 + mas connecting two spins, i.e. 
Q,fc,m = cr = {SqS^) = I (SoSr). Using the NN correlators ci,o,o and co,o,i the GS energy 
per spin Eq/N is given by Eq/N = 3 J||Ci,o,o + 3J^Co,o,i/2. The Eqs. (ITTj) - (|^ contain eight 
parameters Ci,o,o, Ci,i,o, C2,o,0) Co,o,i, Ci,o,i, £0,0,2, Vh ^^^ ^-L' which must be determined self- 
consistently. The relation between the correlators and the corresponding Green's functions 
is given by the spectral theorent^i^. Its application to the correlators Q,^,^ leads to six 
equations, a seventh one is the sum rule Co,o,o = 1/2- One additional equation is obtained 
requiring that the matrix of the static susceptibility Xq^ = Xq^i^ = 0) is isotropic in the 
limit q ^ 0^, i.e. limg^^o,g„^oXq"lq.=o = ^i^g,^oXq'\qa^=qy=o- From that constraint we 



obtain 

Ci,o,o {J± (1 - 2r7_L (3co,o,i - £0,0,2)) + 8J||r7_L (ci,o,i - Co,o,i)) 

-co,o,i (-^11 (1 - 2?7|| (4ci,o,o - ci,i,o - C2,o,o)) +4J± (?7±ci,o,i - ?7||Ci,o,o)) = 0. (23) 

For the discussion of the magnetic GS ordering we consider the spin-spin correlation functions 
(SmaSri/3). A possible magnetic LRO in the system is described via {SmaSnfs) at large 
distances R„ — R^- More precisely, we consider, as in Refs. l44lJ45l . a condensation term 



Cqq/3 in the Fourier transformation Sa/siq) of {Sma^np) according to 

(SmaSn/3) = jr'^Safs (q) exp (-iqr„„,„;3) + -Cq„^ exp (-iQr^a.nfd) , (24) 

where Q is the corresponding wave vector of the magnetic order and Sa/s (q) is given by 

^./5(q) = ^E5^<"l^q)0q|^)- (25) 

The existence of CQa/s 7^ is accompanied by a diverging static susceptibility Xq~ at the 
magnetic wave vector Q. 

To describe the magnetic order in a short-range ordered phase we use in addition to the 



spin-spin correlation functions the correlation length ^^ a 



one a certain direction e^. To 



44lJ45l . We expand that eigenvalue 



calculate C,i, we apply the procedure illustrated in Refs. 
of the static susceptibility Xq ? which is largest at the magnetic wave vector Q around this 
point Q. The relevant magnetic wave vector Q in the short-range ordered phase is that 
vector, where the largest eigenvalue of the static susceptibility has its maximum. Then 
the square of the correlation length along e^, is given by the factor in front of the square 
of the corresponding component of the q- vector q,y = e,y(Q — q). For illustration we give 
the expression for the intralayer correlation length C,\\ in the pure kagome limit {J± = 0), 
^jj = — 2?7||Ci,o,o/ (1 + 27711 (2ci,o,o + ci,i,o + C2,o,o)), and omit the lengthy general expressions 
for the inter- and intralayer correlation lengths. 

V. RESULTS 

A. Antiferromagnetic interlayer coupling J± > 

Solving the set of the eight self-consistency equations we find that there is no solution 
with a non-zero condensation term. Therefore we conclude that within the applied rotation- 
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FIG. 2: First, second and third neighbor intralayer [R = (1,0,0); (1,1,0); (2,0,0)] and interlayer 
[R = (0,0,1); (1,0,1); (0,0,2)] correlation functions in dependence on the antiferromagnetic 
interlayer coupling J_\_. 

invariant Green's function decoupling scheme there is no indication of magnetic LRO for 
the stacked 3D kagome HAFM with antiferromagnetic interlayer coupling. This remarkable 
result is in accordance with the findings within the LSWT (see Sec. IIII|) . To study the 
influence of the interlayer coupling on the GS magnetic short range order we calculate the 
spin-spin correlation functions and the correlation lengths. For convenience we set Jy = 1. In 
Fig. 121 we show several first, second and third neighbor correlation functions. The strongest 
change is seen in the strengths of the interlayer correlators which increase in the interval 
< Jx < 1, for J± K, 1 their values are already close to those of the pure linear chain. 
The change in the intralayer correlators is small. While the NN intralayer correlator shows 
some decrease in strength, the second and third neighbor correlators are weakly increased 
for small J^. Remarkably, already for J^ < 1 the interlayer correlators become stronger 
than the intralayer correlators. Note that in the limit J_l = our values for the GS energy 
as well as for the correlators coincide with available results reported in Refs. |22I J! 



As illustrated above, to calculate the inter- and intralayer correlation lengths ^_l and ^y 
we have to expand the extreme eigenvalue of the static susceptibility around Q, which is 
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FIG. 3: Intralayer (^m) and interlayer (^±) correlation lengths in dependence on the antiferromag- 
netic interlayer coupling J±. 

Q = (0, 0, vr) for all J± > 0. The results for the intra- and interlayer correlation lengths 
are shown in Fig. |H1 The intralayer correlation length is of the order of the NN separation 
d = 1 in the pure 2D limit. It increases up to a maximal value of about ,^|j ^ 2.5d for 
J_L ~ 1 and goes to zero for J± ^ oo. The increase of the interlayer correlation length 
is stronger. Though we expect a diverging ^j_ for the isolated chain, it remains finite for 
any J^ > 0. This behavior is a drawback of the Green's- function decoupling, which is not 
able to describe the critical GS with a power-law decay of the correlation functions of the 
linear-chain HAFM. 

B. Ferromagnetic interlayer coupling J± < 

Now we turn to the case J± < which appears to be more complicated from the numerical 
point of view. Starting from J^ = we find a solution of the self-consistency equations with 
zero condensation term until J± ~ —J\\ = —1. That means we can conclude that for 
ferromagnetic interlayer coupling there is no indication for magnetic LRO for \Jj_\ < 1, too. 
Beyond this point our numerical procedure fails because of numerical instabilities in the 
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q space integrals. The reason for that consists in the specific behavior of the dispersion 
relations c^aq from Eqs. (fTSj) . (f^ and (J22I)- Increasing the strength of the ferromagnetic 
interlayer coupling beyond |J_l| > 1 we find that cuiq becomes more and more a fiat zero 
mode in the Qx-Qy plane at g^ = 0, and at the same time ct;2q and c^sq become more and 
more degenerate leading to divergent terms in q integrals. This behavior of the cj^^q is in 
coincidence with our findings within the LSWT, where we have found a fiat zero-mode and 
two degenerate modes in the g^ — g^-plane at gz = (see Sec. IIII|) . We believe, that the 
physical interpretation of this behavior is connected with the classical limit (large quantum 
number s) of the pure 2D kagome HAFM. In the large s limit the results of the LSWT theory 
become reliable. In the 3D stacked kagome lattice with ferromagnetic J^ the spins along 
the chains become stronger coupled with increasing \J_[_\ and can be considered as large s 
(quasi-classical) composite spins. These quasi-classical spins are coupled kagome-like with 
each other, thus we are most likely faced with the classical HAFM on the kagome lattice 
with its huge non-trivial degeneracy of the GS. Though we do not have rigorous statements 
on the zero-temperature correlation functions of the classical kagome HAFM, presumably 
at T = the average of the spin-spin correlation over the set of all ground states does 
not exhibit magnetic LRO.-^^ Hence we may conclude that the zero-temperature spin-spin 
correlation functions within a kagome layer do not exhibit LRO for any ferromagnetic J^. 

However, the subtile interplay between the tendency to form large s (quasi-classical) 
composite spins and the still remaining (wea k) quan tum fiuctuations may lead to order by 
disorder effects. As it was argued in Refs. I(tI]7| ]53I)54I for the pure 2D classical HAFM on the 
kagome lattice at low temperatures T ^ 0, the entropy of fiuctuations may lead to different 
relative Boltzmann weights of the classical ground states this way favoring planar ordered 
ground states. Presumably there is a nematic LRO but an algebraic decay of the spin-spin 
correlation for T -^ 0^»^. Hence, for the stacked system considered in this paper at T = 
and for J± — > — oo, due to remaining quantum fluctuations a nematic order accompanied by 
an algebraic decay in the intralayer spin-spin correlation functions is possible. 

To describe the magnetic GS in more detail we present, similarly to the case J± > 0, the 
spin-spin correlation functions and the correlation lengths in Figs. |2 El Due to the above 
mentioned numerical problems our data are restricted to < | J_l| < J\\ = 1. Several flrst, 
second and third neighbor correlation functions are shown in Fig0] Again the change in the 
intralayer correlators is very small, whereas the strengths of interlayer correlation functions 
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FIG. 4: First, second and third neighbor intralayer [R = (1,0,0); (1,1,0); (2,0,0)] and inter- 
layer [R = (0,0,1); (1,0,1); (0,0,2)] correlation functions in dependence on the strength of the 
ferromagnetic interlayer couphng |J_l|. 
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FIG. 5: Intralayer (^m) and interlayer (^^) correlation lengths in dependence on the strength of the 
ferromagnetic interlayer coupling |J_l|. 
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increase with growing | J^ | and become larger than those of intralayer correlation functions 
except for the NN correlation. 

The correlation lengths for J± < belong to Q = (0, 0, 0). We show the intra- (^y) and 
interlayer (,^_|_) correlation lengths in Fig. El ^\\ and ^_l increase in the parameter region 
shown but are still of the order of the NN separation d. Obviously, the increase in ^_|_ (i.e. 
along the chains) is significantly stronger. We expect that ^_l diverges for large |J_l|. As 
discussed above, also a diverging ^y for J^ ^ — oo would be possible indicating an algebraic 
decay of the correlation functions. 

VI. SUMMARY 

In this paper we have investigated the GS of the stacked three-dimensional kagome spin 
half Heisenberg antiferromagnet for ferromagnetic and antiferromagnetic interlayer coupling. 
To study the magnetic GS order we have applied a second-order Green's function decoupling 
scheme going beyond the RPA. Though one could expect that an increase in the dimension 
from two to three would stabilize magnetic order, we find some evidence that the interlayer 
coupling is not able to create magnetic LRO within the antiferromagnetic kagome layers. 
These findings based on the Green's function scheme are supported by linear spin wave 
theory. 
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